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Abstract 

The neutron skin and halo distributions in medium and heavy nuclei are calculated within the 
extended Thomas-Fermi approximation. Calculations are carried out for the effective Skyrme-like 
forces using the direct variational method. The analytical expression for the isovector shift of the 
rms radii Ar„p as a sum of skin- and halo-like terms is obtained. The contribution of halo and skin 
terms to Ar„p are found to be approximately equal. 
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I. INTRODUCTION 


For the analysis of the experimental data the simplest two-parameter Fermi (2pF) distri- 
bntion is often nsed: 


F,{r) = 


1 , 

1 -|- exp ' 


-1 


( 1 ) 


where Rq is the half-density radins and Og is the diffnseness parameter of the distribntion. 
Here g = n is for nentron and q = p for proton distribntions. There are two opposite 
pictnres in description of a two-component hnite Fermi system with 2pF. The hrst one is 
the so-called ’’nentron skin-type” distribntion having the nentron half-density radius larger 
than the proton half-density radius, Rn > Rp, and equal diffuseness parameters an = ap. The 
second one is the ’’neutron halo-type” distribution having = Rp and an > ap. A mixture 
of the neutron ” skin-type” and ” halo-type” distributions having Rn > Rp and > ap is 
also possible. It is found that the experimental data can be reproduced by a variety of these 
2pF distributions with different values of — Rp and an — ap [^ . In that sense there is an 
ambiguity in theoretical description of the experimental data. 

In the present paper we study this ambiguity within the extended Thomas-Fermi approx¬ 
imation (ETFA) based on the direct variational method. The nucleon densities Pq{r) are 
generated by the prohle functions which are obtained from the requirement that the energy 
of the nucleus should be stationary with respect to variations of these prohles. 

II. THE MODEL 


In gener^ the total energy of a nucleus includes the kinetic and potential energies and 
is given by 2H4 

Etot = J [ekin(r) + epot(r)] , (2) 

where the potential energy density, epot(r), includes the NN-interaction, eNN(i'), the spin- 
orbit part of the NN-interaction, eso(r )5 and the Coulomb energy density, ec(r): 


epot(r) = eNN(r) + eso(r) + ec(r)- 


( 3 ) 


In the framework of ETJFA, the kinetic energy density is given by the sum of the neutron 
and proton contributions 


ekin(r) = ekin,n(r) + ekin,p(r), 


( 4 ) 
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where 


ekin,qirj = 


2m 


+ -4—+ 

5 ^ 6 b Pq 6 


(5) 


Here m is the bare nucleon mass. The semiclassical consideration gives the value of the 
parameter = 1 in Eq. ([5]) j^, 3|. In the asymptotic limit r —)■ oo, the semiclassical particle 
density Pq with p = 1 falls off to zero significantly faster than the one from the quantum- 
mechanical calculation, where one has p = A. We will use both values of the parameter p to 
study the neutron halo and skin appearances and their superposition in nucleh 
For the effective nucleon-nucleon interaction we will use the Skyrme force 
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where p = -|- Pp is the total density of nucleons. The spin-orbit part of the Skyrme force 

is written as 

m* (r) ITq^ 


esoirj = 


[Pn Pn + Vpp) -|- Pp (Vp„ -|- ‘2'S/Pp) ] . 


(7) 


t)? 4 

Here to, H, t2, ts, xq, xi, X2, X3, z/, Ho are the parameters of the Skyrme interaction and the 
effective mass of a nucleon is 


m*{r) = 


m 


^ , where f3 = 

1 -I- /5p(r) h? 4 


--(3ti -|- 5t2) + t2X2 
4 


The Coulomb energy density is taken in the well-known form { 4 ] 


ec(r) = eVp(r)- / dt 


/ Pp(r') 


-e I - 

TT 


1/3 


Pl’(r). 


( 8 ) 


The unknown values of Pg can be evaluated from the condition of equilibrium. The 
corresponding condition implies that the total energy Etot reaches a minimum value for a 
given number of neutrons, N, and protons, Z, 


N= drpn{r), Z= drpp{r). 


(9) 
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Following the direct variational method the trial function for Pq(r) is taken as a power 
of the Fermi function 


Pq{r) = PO,, 


where po,g 5 Rqi Og and ^ are the variational parameters. 


1 , 

1 + exp ' 




( 10 ) 


III. NUMERICAL RESULTS 


The nuclear rms radius for nucleons is dehned as 



dr pq{r). 


( 11 ) 


The difference between the neutron and proton root mean square radii gives the neutron-skin 
thickness 




( 12 ) 


From the dehnition of trial functions it is seen that the value of Ar^p can be described by 
the different radii Rq (skin effect) and the different diffusenesses (halo effect) of neutron 
and proton distributions, see also jo]. Within the leptodermous approximation aq/Rq -C 1 
one can separate the above contributions and the value of Ar„p is written as 


Ar„p = Arf- + + O , (13) 

where Ar*p*"' and Ar'^p° are caused by the skin and halo effects, respectively. The corre¬ 
sponding values are given by (see Appendix A) 

+ l ['‘o(0 - 2 ki(0] (I)' 

- 5 (78«:i)(0 - 204«:„(f)Ki(?) + SlK^)?)) (|)d , (14) 

and 

+ 1 (75Ki>(?) - 204 k„(Oki(?) + 81k 2(0) (|)d . (16) 

where we have introduced the notations: R = {Rn + Rp)/2, a = {an + ap)/2, — Rp, 

and Aa = an — ap, and /lj( 0 the generalized Fermi integrals (see Appendix A). 
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Fig. 1. The partial contribution to the isotopic shift of radii Ar„p from the skin effect, Ar®p™/Ar„p 
(solid lines) and the halo effect, ArJ(p^°/Ar„p (dashed lines) versus the parameter t]. The calcula¬ 
tions have been performed with SkM parametrization for nuclei ^^*^Sn and ^^®Pb. 

In Fig. [T] the partial contributions to the neutron-skin thickness from the skin effect, 
Ar®p°/Ar„p, and the halo effect, Ar^p°/Avnp, are shown. One can see from Fig. [1] that 
the halo partial contribution increases and the skin partial contribution decreases with the 
increase of parameter r]. Close to the semiclassical approach region, 1 < r] < 2, the partial 
contributions from the halo effect are negative and therefore the partial contributions from 
the skin effect are greater than the unity. The reason for that is the high value of the 
parameter ^ for the asymmetric Fermi function, (see below Fig. [2]), so the densities fall 
off too quickly in the outer surface. This means that the halo component is negative and 
reduces the difference between the neutron and proton root mean square radii. In the 
region 1 < rj < 2 one has more preferably proton halo effect as reported in |^. For the 
values of ?7 ~ 4 — 5, where the tails of the densities are close to that obtained in the 
quantum mechanical calculations, the halo and skin components of Ar„p are found to be 
approximately equal. 

The sensitivity of the partial contributions Ar®p™/Ar„p and Ar^‘^°/Avnp on the Skyrme 
force parametrization is shown in Fig. |3l It is seen from Fig. |3]that the calculated values 
change of about 20% for different parametrizations. 

The dependence of Ar„p on the asymmetry parameter X together with the experimental 
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Fig. 2. Parameter ^ versus r/. The calculations have been performed for nuclei ^^^Sn and 208 p]^ 
using SkM parametrization. 


data are shown in Fig. 01 As it is seen from the hgure the calculation with the ’’quantum 
mechanical” value rj = 4 (solid curve), in general, agrees better with the experimental data 
than for t] = 2 (dashed curve). The obtained dependencies are approximately linear for the 
considered values of isotopic asymmetry parameter X. 

In Fig. [5] the analogous calculations for tin isotopes together with the experimental data 
are shown as a function of mass number. Similarly to the previous cases, the calculation 
with 1 ] = 4 describes the experimental data better than the one with rj = 2. 

The neutron-to-proton density ratio in the peripheral region with respect to the bulk 
value is described by the peripheral halo factor js]. Following j?, ^ we will also use the 
theoretical estimation of the halo factor as 


/■halo / N ^ Pnjr) Z 

ItheoAr) ^ (lb) 

In Fig. |6] we show the calculations of the theoretical halo factor for two values of the 
parameter rj together with the experimental data [^. As in the previous hgures the solid 
lines correspond to the calculation with r] = 4 and dashed lines correspond to rj = 2. The 
left panel is the case for tin isotope ^^^Sn, and the right panel is for ^°®Pb. It is seen from 
Fig. [6] that the calculation with rj = 4 describes the experimental data better than the one 
with 7] = 2. 
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Fig. 3. The partial contribution to the isotopic shift of radii Ar„p from the skin effect, Ar®p™/Ar„p 
(solid lines) and the halo effect, ArJ(p^°/Ar„p (dashed lines) versus the parameter t]. The calcula¬ 
tions have been performed for ^^*^Sn with the following Skyrme force parametrizations: 1 - Sill, 2 
- SkM, 3 - SLy230b. 



Fig. 4. Difference Ar„p between the rms radii of the neutron and proton distributions as a function 
of A = (N — Z)IA. The solid and dashed lines are calculations with r? = 4 and 2 for SkM 

Q 

parametrization of the Skyrme force. The experimental data are taken from Ref. [7|]. 
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Fig. 5. Difference Ar„p between the rms radii of the neutron and proton distributions as a function 


of A for Sn isotopes. The solid and dashed 
The experimental data are taken from Ref. 


ines are calculations with t] = 4 and 2 for SkM force. 


fl. 



Fig. 6. The halo factor in ^^^Sn and ^*^®Pb as a function of the distance from the center of the 
nucleus calculated using the Eq. (fTbll . The values of the halo factor deduced from experiment are 

Q 

marked by crosses [7|]. 

IV. CONCLUSIONS 

We have considered the neutron skin and halo effects for medium and heavy nuclei within 
the extended Thomas-Fermi approximation using effective Skyrme-like forces and the direct 
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variational method. The nucleon densities Pp{r) and Pn(i‘) are generated by the prohle 
functions which are obtained by the requirement that the energy of the nucleus should be 
stationary with respect to variations of these prohles. 

Using the leptodermous properties of the prohle nucleon densities Pp{r) and Pn(r), we 
have obtained the analytical expression for the isovector shift of the rms radii Ar^p as a 
superposition of the two terms. The hrst one, describes the neutron ”skin-type” 

distribution and the second one, Ar^°, describes the neutron ’’halo-type” distribution. 

Numerical calculations show that the partial contributions to the neutron-skin thickness 
from the skin effect, Ar®p°/Ar„p, and the halo effect, ArJjp°/Ar„p, depend on the parameter 
p. With the inscrease of p from 2 io p = A the partial contributions to the neutron-skin 
thickness from the skin effect decreases from 100% to about of 50% and the halo component 
increses from zero to about of 50%. 

The calculated isovector shift of the rms radii Ar„p and halo factor /theor(^) satis¬ 
factory description for the corresponding experimental data at ?7 = 4. In this case we have 
superposition of the neutron skin and halo effects with approximately equal contributions. 
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APPENDIX A 

In this Appendix, we will consider the rms radius for nucleons as 




f drr^pg(r) 

fdr pg(r) 



(Al) 


where 



Whithin the leptodermous approximation Qg/Rq << 1 one has 


(A2) 


p3 / / \ 2 / \ 

l + 3/.o(0^ + 6«:i(0(^) +3«^2(0(^) , (A3) 

O \ ihq y itq J y ihq J j 

p5/ f /\^ f 

1 + 5«:o(0^ + 20«:i(o(^) +30«:2(o(^) +20k,{o{^) 

O \ Ihq y Ihq J y ihq J y ihq J 

+ 5«:4(0f^y), (A4) 

y / j 
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where Kj{^) are the generalized Fermi integrals derived in Ref. 

POO 

«,( 0 = / dxx^[{l + en-^-{-iy + 

Jo 

Inserting Eqs. flASp and flA4D into the Eq. flAip one gets 

- vfis, {i + ^ - 2 ki(0) (^) 

\ (75ko(0 - 204ko(0«^i(0 + 81^2(0) 
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From Eq. flA6[) the isovector shift of the rms radii Ar„p reads 

^rnp = VW) - \fW) 

~ '\/f -Rp + «o(0(an -ap)-^ - 2ki(0) 

+ g — 204 ko( 0 ^i (0 + 81^2(0) 

One can rewrite last terms of Eq. flA7p using the relations 


<_4 


(|4^) + 2 

\ “r / 


dfL / \ 

Dip ~ 

■P^n ~r -Ptry 


_ yp_ ^ _2 ( 

R2 i?2 - 
•'n V 


Rn ~l~ 


(K - Rp) + 3 ^ ") (a„ - ap) 

\ “r -thn j 


Finally, we have 


Ar„p - V 5 11 + 2 (rI + 4 


(AS) 


(A6) 


(A7) 


(AS) 

(A9) 


- ^ (75Ki)(0 - 204«:o(0«^i(0 + 81^2(0) () } (Rn - Rp) 

o y -t^n “r -thp ^ 


(AlO) 


+ |^o(0 “ 7 (ko(0 ~ 2ki(^)) 


dfi Qjp 

Rn H“ Rv 


^ (75ko(0 - 204ko(0«;i( 0 + 81fi:2(0) Q"" | (on - ap) 


(All) 
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